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Abstract. Off-diagonal profiles </>cid(^) of local densities (e.g. order parameter 
or energy density) are calculated at the bulk critical point, by conformal methods, 
on a strip with transverse coordinate v, for different types of boundary conditions 
(free, fixed and mixed). Such profiles, which are defined by the non-vanishing matrix 
element {0\(f>{v)\{f>) of the appropriate operator (p{v) between the ground state and the 
corresponding lowest excited state of the strip Hamiltonian, enter into the expression 
of two-point correlation functions on a strip. They are of interest in the finite-size 
scaling study of bulk and surface critical behaviour since they allow the elimination 
of regular contributions. The conformal profiles, which are obtained through a 
conformal transformation of the correlation functions from the half-plane to the strip, 
are in agreement with the results of a direct calculation, for the energy density of the 
two-dimensional Ising model. 



Following the pioneering work of Fisher and de Gennes [1], the study of order 
parameter and energy density profiles near surfaces has been an active field of research 
during the past years. These profiles have been calculated at, and near the critical 
point, in the mean-field approximation [2], using field-theoretical approaches [3] 
and through exact solutions [4, 5]. Much progress has also been achieved in their 
calculation at bulk criticality in two-dimensional (2d) systems making use of conformal 
techniques [6-12]. 

In a semi-infinite 2c? system, the profile 4>{y) of a fluctuating quantity, such as the 
order parameter or the energy density, is obtained in the transfer matrix formalism as 
the diagonal matrix element (O|0(?;)|O) of the corresponding operator (j) in the ground 
state |0) of the Hamiltonian H = — InT, where T denotes the transfer operator along 
the surface. 

On a strip of infinite length and finite width L, one may also consider the off- 
diagonal profile 4>^^{v), where v is the transverse coordinate. The profile is then defined 
as the off-diagonal matrix element (O|0(u)|(/)) between the ground state |0) of the 
Hamiltonian Ti. on the strip and its lowest excited state |0) leading to a nonvanishing 
matrix element. 
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These off-diagonal profiles are commonly used at the bulk critical point to obtain 
informations about the siu'face and bulk critical behaviour via finite-size scaling, while 
avoiding the regular terms which contribute to the diagonal ones. Actually, in the 
absence of an external symmetry-breaking field, an off-diagonal matrix element has 
to be used to study the scaling behaviour of the order parameter since the diagonal 
one vanishes, due to symmetry. They are also of interest because of their high degree 
of universality. One may mention two recent studies where off-diagonal profiles are 
considered at bulk criticality: the spin 1/2-spin 1 Ising quantum chain in [13] and the 
random Ising quantum chain in [14]. 

In the following, the scaling form of such off-diagonal matrix elements is obtained 
at the bulk critical point for two-dimensional conformally invariant systems. It 
is deduced from the asymptotic behaviour of the appropriate two-point correlation 
function in the half-space after a conformal transformation to the strip geometry. We 
consider first symmetric and then asymmetric boundary conditions on the strip and 
compare our results to some exact expressions of the energy-density profiles, obtained 
in the appendix, for the Ising model in the extreme anisotropic limit. Since we use 
conformal methods noninvariant boundary conditions, like a finite surface field for 
which an interesting scaling behaviour has been recently observed [15, 16], and off- 
critical systems are here excluded. 

Let us first briefiy review the Fisher-de Gennes result and its consequence for 
the profile on a strip. We consider a nonvanishing profile 4>{y) on a semi-infinite 2d 
conformally invariant sytem at its critical point with a surface at y = 0. It may be 
the energy density profile with any type of uniform boundary conditions or the order 
parameter profile with fixed boundary conditions. The problem involves a single length 
scale, the distance y from the surface. Under a length rescaling by a factor b, the profile 
transforms as (f){y/b) = b'^'i'(f){y) so that 

<t>{y) = Ay-^^ (1) 

where is the bulk scaling dimension of (p. Now, making use of the conformal 
transformation w = {L/Tr)\nz, with z = x + iy = pe'^ and w = u + iv one obtains 

p = exp(-), 0=-, (2) 

and the local dilatation factor is b{z) = \dw/dz\~^ = wp/L. The half-plane y > 

is transformed into a strip — oo < u < +oo, < ?; < i with the same boundary 
conditions as the half-space on both edges [17]. The profile transforms as [18]: 

(l>{w) = b{zr*^{z) (3) 

so that [6] 



L . /"KV 

— sin I 

TT V 



Trv\ 

t) 



(4) 



One may notice that the surface critical behaviour is hidden in the sine variation. 
At a fixed distance I <C L from the surface the profile behaves as 



<f>{l)^Al-^^{l + ^n^l'x^L-^ + ---) , 



(5) 
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in agreement with the Fisher-de Gennes conjecture which gives a L correction 
term [1] with an universal ampUtude [9]. Generally, the exponent d is expected to 
give the scaling dimension of (j) at any surface transition leading to a nonvanishing 
profile in the semi-infinite critical system [7,9,19-21]. The 0{N) model at the special 
transition (for TV < 1 in 2d.) provides a counter-example, the surface energy exponent 
being smaller than d in this case [22, 12, 23] ||. 

Let us now turn to the calculation of the off-diagonal matrix element. On the 
critical semi-infinite system, conformal invariance strongly constrains the form of 
the connected two-point correlation function, Q^°f^{xi — ajg, 1/1,2/2)- Applying an 
infinitesimal special conformal transformation which preserves the surface geometry, 
one obtains a system of partial differential equations for Q^'^, from which the following 
scaling form is deduced [24]: 

- ^2,yi,2/2) = {yiv^r^'gi^} , ^ = " + " . (6) 

i/li/2 

When w 1, ordinary scaling leads to g{u)) ^ u>~'''i> where a;^ is the scaling dimension 

of <j) at the surface. This limit also corresponds to where, using polar 

coordinates, (6) may be rewritten as: 

G;T{Pi,P2,0„e,) ~ (P1P2)-"* (^) ' (sin ^1 sin ^2)"*""* • (7) 

We now consider the connected two-point function at bulk criticality in the strip 
geometry. In the same limit, it can be obtained through the conformal transformation 
of (7) from the half-plane to the strip geometry with [17] 

^^4^1' ^2) = K^l)"*K^2)"*^^^(^l>^2) . (8) 

giving: 

GlTi^i - u„v„v,) = [^y^* {PiP2r'G;TiPi^P2,0^,9^) 



/ TT \ 2x^ TTX^ . I . /TrV^\ / TTV2 \ 



'-Ip -^v 



On the strip, making use of the transfer operator e~^, it can also be written as an 
expansion over the eigenstates \n) of the critical Hamiltonian H. with eigenvalues 
and ground state energy E^: 

GlTiui - U2,v„v^) = ^{0\4>M\n){n\^{v^)\0)exp[-{E^ - Eo){u^ - u^)] . (10) 

n>0 

In the limit ^ P2 which corresponds to Ui ^ U2 on the strip, the sum is dominated 
by the contribution of the lowest exited state ji^) with a nonvanishing matrix element 
so that: 

GlTi^i - U2,v„V2) ^ {0\4>{v^M{mv2meM-iE^ - E^){u^ - u^)] • (11) 



II We thank E Eisenriegler for pointing out this exception to us. 
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Comparing with (9) we recover the gap-exponent relation E, — Eq 



l/L [17] as a 



by-product and we may identiiy the expression of the off-diagonal matrix element 



sm 



/TTV 



(12) 



with symmetric boundary conditions on the strip. The off-diagonal energy-density 
profile for the Ising model given in equation (A5) of the appendix, which is valid for 
symmetric free or fixed boundary conditions, agrees with (12) in the continuum limit. 

At a fixed distance I <C L from the surface, 4>odi^)' behaving as L~^'i>, is quite 
appropriate to perform a finite-size scaling study since the regular term which appeared 
in (5) is now avoided. One may notice that the diagonal profile (^{v) in (4) is formally 
recovered with = in (12). Finally, in the half-plane limit (L ^ oo), the amplitude 
of the off-diagonal profile vanishes. This is consistent with the scaling considerations 
leading to equation (1): a nonvanishing profile on the semi-infinite system necessarily 
gives the diagonal profile (4) on the strip. 

Next we consider the case of mixed boundary conditions on the half-space. This 
means having different scale-invariant boundary conditions a and b on the positive 
and negative x-axes, respectively, with a,b — f (free), -|- or — (fixed). The plane- 
to-strip conformal transformation leads to different boundary conditions a and b 
on opposite edges of the strip. The method employed to determine the scaling 
behaviour of the two-point function in (6) can no longer be used since it supposed 
translation invariance along the x-direction which is now broken by the mixed 
boundary conditions on the surface. Fortunately, the two-point correlation functions 
with mixed boundary conditions have been obtained explicitly for the Ising model 
using conformal methods [11]: n-point correlation functions in the half-space with 
mixed boundary conditions arc determined by the same differential equation as a 
particular 2n + 2-point bulk correlation function. Using the results of Burkhardt and 
Xue [11], when Pi/p2 ^ 1, the asymptotic spin-spin correlation functions read 



(PiP2sin6'iSin6l2) 



cos 01 COS 02 + sin^ 6-^ sin^ ^2 ( ~ ) + 



Oa^ ~ (Pl/02Sin^l sin e2)"'/'(cOS ^1008^2)'/' 



1 -|- - sin 9-1 sin ^2 tan tan ^- 



2 \ P2 



(13) 



2; Pi 



+ • 



whereas the following expressions are obtained for the energy-energy correlation 
functions: 



(P1P2 sin^li sin6l2)" 



St/ ~ (PiP2sin6'isin6l2) ^ 



(l-4sin2 6li)(l-4sin^ 6*2)+ 



+64 sin^ 61 sin^ 62 cos 9^ cos ^2 



cos 9i cos ^2 + 8 sin^ 9^ sin^ ^2 ( ~ 



(14) 
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Making use of the conformal transformation to the strip geometry as in the first hne 
of c;quation (9) and c;omparing to the expansion of the correlation function on the 
strip with asymmetric boundary conditions, one may identify the profiles. Since the 
correlation functions in (13) and (14) are the unconnected ones, the expansion in (10) 
now contains the term n = 0. In the limit u-^ — ^2 3> 1 the leading contribution is the 
product of the diagonal profiles 0(t;^)0(t;2) and the next term contains the product of 
the off-diagonal ones, as before. In this way, one obtains 




^:»-^(f)"'*(^). 

for the off-diagonal order parameter and energy density profiles. The results of a direct 
calculation of the off-diagonal energy density profiles, in equations (A6) and (A9) of 
the appendix, are in agreement with the two last equations of (15) in the continuum 
limit. 

Diagonal profiles on strips with asymmetric boundary conditions for the Ising, 
Q-state Potts and 0{N) models can be found in references [8,10-12] where they were 
deduced from the appropriate one-point functions in the half-plane. 

With mixed boundary conditions the correction term in p^j gives a first gap 
vanishing as tt/L on the strip: the connection with surface exponents we had before is 
lost. Contrary to the case of symmetric boundary conditions in equations (4) and (12), 
there is no simple functional form adapted to the different boundary conditions. The 
diagonal profiles generally contain a regular leading contribution in their surface 
finite-size scaling behaviour whereas the off-diagonal ones behave as L~(^*)<»''' with 
(x^)+,_ = (4)+^. = (x|)^. = 2 and [xDj - 1/2. 

One may notice that diagonal and off-diagonal profiles obtained in [14] for the 
random Ising quantum chain are in quite good agreement with the conformal results 
given above although the system, which is strongly anisotropic, is not conformally 
invariant. It would be interesting to check wether this is peculiar to the Ising model 
or if it holds true for other Id anisotropic systems as well. 



Appendix. 

In this appendix, we calculate the off-diagonal energy density profiles on a strip for the 
Ising model with different boundary conditions. We work in the extreme anisotropic 
limit [25] where the critical Hamiltonian reads 



n 



L-1 



L-1 



J=2 



(Al) 
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where <J^{j) and <J^{j) are Pauli matrices at site j. The diagonization proceeds in 
two stops [26, 27]: first the Hamiltonian is rewritten as a quadratic form in fermion 
operators, c^(j) and c(j), using the Jordan- Wigner transformation [28] and then it is 
diagonaHzed through a canonical transformation to new fermion operators, jjI and r;^ , 
such that: 



■ M)Wk + ifkU) - Mj)]vk 



(A2) 



iff. and tpi^ arc normalized eigenvectors, which are conveniently calculated by a 
diagonalization process, as described in [29] . 
With free boundary conditions, h^= hr 



fkij) = (-1)^ 



%/2L + l 
2 



cos 



L - 

2k 



1, the eigenvectors are given by 

1 /j-r 



2L+1 

2k -1 



: sm 



with fc = 1, 2, 



y/2L+l "'\2L + V 
. , L. The off-diagonal energy density profile is expressed as 



{j) = {0\a^U)\e) = 2{0\c^j)c{j)\e}, 



(A3) 



(A4) 



where the ground state ]0) is the fermion vacuum and \e) is the lowest eigenstate with 
two fermionic excitations so that one obtains: 



'od 



(i) 



4 



2L+1 



sm 



2L-I-1 



cos 



^^i+ll^] -cos( ^ \m\^^t:^]V^^^ 
2L+1 \ ''°\2L+lj'' [ 2L+1 Jj- 

When L ^ 1, the first term may be neglected and the second one gives the off-diagonal 
profile in equation (12) since, for the 2d Ising model, = 1 and a;^ = 2 at the ordinary 
surface transition. 

When the two boundary spins are fixed one has to put = hj^ = and then the 
Hamiltonian in equation (Al) describes both the -|— |- and the H — boundary conditions. 
Instead of solving this problem explicitly one may use the duality properties of the 
quantum Ising model [25] . Through a duality transformation, the critical system with 
-|— I- (-1 — ) boundary condition is related to the energy (magnetization) sector of the 
critical free chain^, with the correspondence L — > L — 1 and j — > j — 1. Then the 
^od U) energy density profile is directly given by (A5) with the above substitutions. 

To calculate the -I — profile one considers the magnetization sector of the free 
chain, where the two lowest states are r?|]0) and tjIIO). Then the ^^^"(i) profile is 
given by: 

4d"0') = v>iU)M) + v>2U)M) 

16 



2L-1 



■ cos 



r(2j"l)7r1 


sin 


r(j- 1)7^1 


•|cos^ 


r(2i-l)7rl 


— sin^ 




} 


4L-2 


2L-1 


4L-2 


2L-1 





1 The energy (magnetization) sector corresponds to an even (odd) number of fermion excitations in 
the system. 
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which asymptotically behaves as £^^{j) ~ 2sm{2jn/L)/L in agreement with the 
conformal result in (15). 

The +/ boundary condition is realized with a vanishing transverse field = 
on the first spin. Then ^^.(l) commutes with H and remains fixed in one of its 
eigcnstatcs corresponding to crj.(l) = ±1. This introduces a vanishing excitation in 
the system which doubles the spectrum. The behaviour of the zero-mode eigenvectors 
is anomalous: 



Vo(i) = o, ^o(i) = (-i)^+y^- 

The next fermion state is characterized by the eigenvectors: 

(i - l)7r ' 
L 

(j - l/2)7r' 



(A7) 



V'i(i) = (-iy+S/;^cos 



(A8) 



In terms of eigenvectors the profile is similar to (A5) and, due to the vanishing of 
(^q(j), it takes the form 



4d (i) = -'PiiMoU) = ^sin 



(A9) 



in agreement with (15) in the continuum limit. 

The results in equations (A5), (A6) and (A9) have been checked numerically on 
small-size quantum Ising chains. 

We thank H Riegcr, M Henkel and D Karcvski for helpful discussions and 
E Eisenriegler for a critical reading of the manuscript. FI is indebted to the Hungarian 
National Research Fund for financial support under grants no OTKA TO12830, 17485 
and 23642. 
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